Local Magnetic Inhomogeneities in Lightly Doped BaFe 2 As 2 



A. P. Dioguardi, 1 N. apRoberts- Warren, 1 A. C. Shockley, 1 S. L. Bud'ko, 2 N. Ni, 2 P. C. Canficld, 2 and N. J. Curro^E 

1 Department of Physics, University of California, Davis, CA 95616 
2 Ames Laboratory U.S. DOE and Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011, USA 

(Dated: August 19, 2010) 

We report 75 As NMR measurements in BaFe2As2 doped with Ni. Like Co, Ni doping suppresses 
the antiferromagnetic and structural phase transitions and gives rise to superconductivity for suf- 
ficiently large Ni doping. The spin lattice relaxation rate diverges at Tn, with a critical exponent 
consistent with 3D ordering of local moments. In the ordered state the spectra quickly broaden 
inhomogeneously with doping. We extract the average size of the ordered moment as a function of 
doping, and show that a model in which the order remains commensurate but with local amplitude 
variations in the vicinity of the dopant fully explains our observations. 
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The recent discovery of the iron arsenide superconduc- 
tors has reignited interest in the interplay of supercon- 
ductivity and magnetism in condensed matter*^ There 
are several iron arsenide families that exhibit supercon- 
ductivity cither under pressure or with chemical doping. 3 
Each family, however, contains a common structural el- 
ement consisting of FeAs planes, in which the Fe 3d or- 
bitals hybridize with the As p orbitals, giving rise to mul- 
tiple bands that cross the Fermi energy^£ In the parent 
state (either undoped or at ambient pressure) the nesting 
of two of these Fermi surfaces gives rise to a spin den- 
sity wave instability^ It is believed that doping (both in- 
plane and out-of plane) and/or pressure tunes the chem- 
ical potential and modifies the nesting condition, which 
consequently suppresses the long range antiferromagnetic 
order. For sufficiently large electron or hole doping su- 
perconductivity emerges.— Although it appears that the 
superconductivity may be intimately related to the anti- 
ferromagnetic fluctuations present in these systems, de- 
tails of how these fluctuations emerge and are related to 
the superconductivity remain clouded.— Furthermore, it 
is unclear how the superconducting condensate can re- 
main intact in the presence of a high concentration of in- 
plane impurities^ Detailed experimental studies of the 
local effects of dopants therefore are crucial to under- 
stand the physics that drive the evolution of long range 
order in these materials. 

In order to shed light on the influence of dopants on 
the antiferromagnetic order, we have conducted 75 As Nu- 
clear Magnetic Resonance (NMR) studies in a series of Ni 
doped Ba(Fei_ a; Ni x )2As2 crystals. The advantage of this 
particular system is that large high quality single crystals 
are available, and Ni has the greatest effect of any transi- 
tion metal dopants so that superconductivity is reached 
at only ~2%^1£ We find that the As NMR spectrum in 
the antiferromagnetic phase broadens inhomogeneously 
very quickly with doping, reflecting a large distribution 
of local hyperfine fields. This distribution can be under- 
stood by realizing that the local spin polarization sur- 
rounding the dopant site is not only reduced on-site (at 
the Ni), but also extends to several of the surrounding 
Fe sites. In the paramagnetic state, we find that the 




FIG. 1. (color online) The local structure of the hyperfine 
fields at the As surrounding a Ni impurity. Black (blue) ar- 
rows indicate the direction of the ordered moment (As hy- 
perfine field). The hyperfine field at the nearest neighbor As 
sites acquires components in the ab plane. 

nuclear spin-lattice-rclaxation rate reflects a continuous 
evolution of three dimensional critical spin fluctuations 
associated with the phase transition at Tn(x). 

Single crystals of Ba(Fci_. T Nia;)2As2 were grown by the 
self-flux method as described inZ*^. The Ni concentra- 
tions were determined via microprobc analysis j 7 -^ NMR 
spectra were obtained by integrating the spin echo as a 
function of applied external field oriented along the crys- 
tal c direction at constant frequency / = 48.28 MHz. 
There is a single As site per unit cell located symmetri- 
cally between four nearest neighbor Fe atoms, with stag- 
gered displacements along the c-axis (see Fig. [TJ. 75 As 
has spin I = 3/2, and the resonant condition is given by 
the nuclear spin Hamiltonian: 



H = -yhI x H 



hv c , 



3Ii 



■n hl , (i) 



where 7 = 0.7292 kHz/G is the gyromagnetic ratio, Hq 
is the external applied field, I a are the nuclear spin oper- 
ators, v cc is the component of the electric field gradient 
(EFG) tensor along the c-direction, 77 is the asymmetry 
parameter of the EFG tensor, and T-Lm is the hyperfine 
interaction between the As nuclear spins and the Fe elec- 
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FIG. 2. (color online) Field-swept spectra of the 75 As in 
Ba(Fei_ a; Ni: C )2As2 at constant frequency of 48.28 MHz for 
various dopings and temperatures. The data for x = 2.6% 
at 20K has been shifted down by 10 K for comparison with 
other dopings. Gray dashed (dotted) lines indicate T 3 (Tn)- 
The solid black lines are simulations as described in the text, 
for q = 0.6 and A = 2.5a. 

tron spins. The hyperfine coupling is given by: 



(2) 



where the sum is over the four nearest neighbor Fe spins 
S(rj), and the components of the hyperfine tensor B are 
given by: B aa = B bb = 6.6 kOe//i B , B cc = 4.7 kOe//i S 
and B ac = 4.3 kOc//ia as determined in the parent com- 
pound by Kitagawa et at— 

Fig. [5] shows field-swept spectra for the As for 
several different dopings as a function of tempera- 
ture. For sufficiently high temperatures in the paramag- 
netic phase Ba(Fei_ x Nia;)2As2 has tetragonal symmetry 
(space group IA/mmm), and hence the EFG at the As 
site has axial symmetry with rj = 0. The spectra consist 
of three resonances at fields Hq = (/ — ?w cc )/7, where 
n = —1,0, or 1, and v cc w 2.5 MHz. v cc is strongly 
temperature dependent^ but surprisingly we find lit- 
tle or no change in the average value of v cc with dop- 
ing. Rather, the spectra reveal an increasing quadrupo- 
lar satellite (n = ±1) linewidth, reflecting an distribution 
of local EFGs as expected in the presence of dopants. By 
x = 3.3% the disorder-induced broadening is so large that 
we are no longer able to discern the quadrupolar satel- 
lites. The increased quadrupolar linewidth, however, is 
minor compared with the large distribution of hyperfine 
fields in the antiferromagnetic state. 

Below a temperature T s (x) the system undergoes a 
structural distortion to an orthorhombic phase (space 
group Fmmm). In the parent compound T s (0) = 134 K 
coincides with TV and the sudden appearance of antifer- 
romagnetic order. Upon doping, however, this first order 
phase transition quickly becomes second order, and the 
two transition temperatures T s (x) > Tn(x) separateiiS 
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FIG. 3. (color online) Hi nt (left axis), Tn and T c (right axis) 
as functions of doping in Ba(Fei_ a; Ni :E )2As2. The open points 
are determined by bulk measurements^ and the solid lines 
are guides to the eye. INSET: The reduced moment M/Mo 
versus reduced temperature T /Tn as measured by NMR in 
Ba(Fei_ a; Ni I )2As2 (□) and by elastic neutron scattering in 
Ba(Fei_ x Co :l; )2As2 (°)-—- The colors correspond to the same 
dopings as in Fig. [5] 

These transitions are reflected in the NMR spectra by a 
change in the EFG at T s (x) as well as the onset of static 
internal hyperfine fields below Tn(x) that shift the spec- 
tra both to higher and lower resonance fields (Fig. [2]). 
In the presence of the internal field Hj n t the resonance 
condition becomes / = 7|Ho + Hi nt | -\-nv cc . Since Hj nt is 
either parallel or antiparallcl to Ho ( 1 1 c) there will be six 
resonances at: Ho = (/ — nis cc )/j±Hi n t, as seen in Fig[2] 
for undoped BaFe 2 As 2 below T N . We find H int = 14.95 
kOe, in agreement with previous work^i This observa- 
tion is consistent with the hyperfine coupling (Eq. [2]) for 
magnetic ordering given by Q = (—,0,0) with moments 
So = 0.87^ B along the (100) direction, as observed by 
neutron scattering (see Fig. [IJJ^ It is crucial to note 
that even though the moments lie along the (100) direc- 
tion, the hyperfine field lies along (001). The reason is 
that the As is located symmetrically between four near- 
est neighbor Fe atoms so that the components of Hi nt in 
the ab plane cancel out by symmetry. A similar effect 
is at play for the planar oxygen site in the cuprates^ 
except that in the iron arsenides the fact that the As is 
not in the plane of the Fe atoms gives rise to the c-axis 
component of the hyperfine field. 

Upon doping, the internal field is reduced and the spec- 
tra broaden dramatically in the ordered state. The spec- 
tra are best characterized by a distribution of H; nt . The 
average internal field (ifint) (determined by the peak of 
the spectra) tracks Tn(x) as a function of doping and 
is shown in Fig. [3J Indeed, H int ~ So, is propor- 
tional to the sublattice magnetization, So, and is a direct 
probe of the antiferromagnetic order parameter. The re- 
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duced internal field, Hi nt (x,T)/ Hi nt (Q,T), is shown in 
the inset. For comparison, recent elastic neutron scatter- 
ing measurements on Co doped Ba(Fei_^Co x )2As2 are 
shown as wellJ^ The similar trends in the doping and re- 
duced temperature matches those of the phase diagrams 
of these two systems.—^ For sufficiently large dopings 
(x > 2.6%) the width of the internal field distribution 
becomes so large that the spectral features are washed 
out. Similar effects were observed in Ba(Fci_ x Co x )2As2, 
and attributed to the presence of incommensurate mag- 
netic orderj 15 ' 16 Although the spectra reveal a broad dis- 
tribution of Hj nt j it is unclear whether there is an inco- 
herent distribution of local Fe spin polarizations, S(r), or 
an incommensurate modulation. Both effects would give 
rise to a distribution of internal fields consistent with 
the spectra. However, in the strongly correlated antifcr- 
romagnctic state a local impurity (Ni) can give rise to 
dramatic changes to the local spin structure even in the 
presence of commensurate order. Kemper and collabora- 
tors have calculated the local spin density in the vicinity 
of a dopant atom in the antifcrromagnetic state of the 
BaFe2As2 system using density function theory (DFT), 
and find that the polarization of the neighboring Fe are 
strongly modified as result of the impurityJ^ Since the 
hyperfine interaction (Eq. [5]) couples each As nucleus to 
four Fe sites, small changes to the spin of any Fe site can 
alter the hyperfine field of several As nuclei. If one of 
the Fe spins were slightly reduced due to its proximity 
to a dopant Ni (or Co) atom, then the hyperfine field at 
several neighboring As sites would acquire a component 
in the ab plane (see Fig. [1} and the resonant fields would 
be given by: 

#o = ^{f-nu cc y-(H^Y HP nt . (3) 

Clearly, a distribution of cither H-^ t or Hf ul could 
broaden the NMR spectra significantly. 

In order to quantify the effect of perturbations of the 
local hyperfine field in the vicinity of a dopant atom, we 
have simulated the spectra using a simple model designed 
to mimic the spin density observed in DFT calculations. 
We start with a 2D square lattice of 100 x 100 Fe spins 
with impurities located randomly at positions r; mp . The 
spin moments are given by S(r) = So$a;(r) cos(Q • r), 
where ^(r) = 1 — aj^exp [— (r — r ; ) 2 /2A 2 ] represents 
the reduction of the local spin polarization in the vicinity 
of a dopant, a is a constant that represents the suppres- 
sion of spin density magnitude at an impurity site, and 
A is a constant that represents the spatial extent of the 
suppression. This model captures the relevant features of 
the DFT calculations, which indicate that not only is the 
impurity site spin moment renormalized, but the magne- 
tization over several neighboring sites in the Fe plane is 
also reduced^ Given this spin density profile, we cal- 
culate the hyperfine fields at each As site in the lattice 
using Eq. [51 and determine the histogram of resonant 
fields using Eq. [3J These histograms are shown in Fig. [5] 
for a — 0.6 and A = 2.5a. The simulated spectra clearly 
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FIG. 4. (color online) (TiT) -1 versus temperature in 
Ba(Fei_ a; Ni a ;)2As2 for x = (■, gray, reproduced froroii), 
x = 0.72% (•, blue), x = 1.6% (▼, green), x = 2.6% (♦, yel- 
low), and x — 3.3% (A, red). Solid lines are fits as described 
in the text. INSET: (TiT)\ F versus T using the same sym- 
bols as in the main figure. Solid lines are linear fits indicating 
behavior characteristic of 3D fluctuations. 

exhibit the same general trend versus doping as our ex- 
perimental spectra, and the values of a and £ are rea- 
sonably close to the DFT calculations, suggesting that 
the origin of the broad NMR spectra is indeed natural 
consequence of doping-induced disorder in the antifcrro- 
magnetic state rather than a change from commensurate 
to incommensurate ordering. 

It is important to note that this non-local effect of 
the Ni dopants is a consequence of the strong correla- 
tions present in the long-range ordered antifcrromagnetic 
state. In the paramagnetic state, NMR experiments show 
little or no change of the spin fluctuations at As sites 
that are nearest neighbors to the dopant atoms.— In this 
case, the magnetic correlations are not sufficiently large 
to induce the non-local effects seen in other strongly cor- 
related systems . 20 ' 21 

In order to investigate the evolution of the spin dynam- 
ics with doping, we measured the spin lattice relaxation 
rate, , at the central transition (I z = —h ■<-> |). The 
data are well fit to a single T\ component, and the tem- 
perature and doping dependence are shown in Fig. [4] 
(TiT) -1 probes the dynamical spin susceptibility of the 
Fe spins, and is given by: 

— =^k B lim VF|(q) Xg(q,fa,) , (4) 

q,/3 

where the form factor ^(q) is the Fourier transform of 
the hyperfine coupling and x|s( < l> w ) i s the dynamical 
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susceptibility. We find that (TiT) -1 exhibits a dra- 
matic enhancement in the paramagnetic state just above 
Tjv due to critical slowing down of the spin fluctuations.— 
For x = 0, (XiT) -1 exhibits a discontinuity at Tjy, con- 
sistent with the first order nature of the phase transi- 
tion. Under doping, however, the transition becomes 
second order<21 The data are well fit to (TiT) -1 = 
a + b/(T — TV) 1 / 2 , where the first term represents a 
Korringa relaxation by itinerant quasiparticles, and the 
second term arises from 3D fluctuations of nearly anti- 
ferromagnetic local moments^ TV is plotted in Fig. [31 
and agrees well with bulk measurements The inset of 
Fig.H]shows(T 1 T)i F = (l/(r 1 r) 1 - a )- 2 ~ where 
£af is the antiferromagnetic correlation length. Clearly 
Iaf ~ (T - T N )~ y with v = 1/2, the mean field result 
for 3D fluctuations. Upon doping, these critical fluctua- 
tions are suppressed as Tjv(x) decreases. Similar trends 
were observed in Ba(Fei_ a; Co :c )2As2, and attributed to 
the evolution of both inter- and intra-band scattering^ 
For sufficiently large electron doping, the hole Fermi sur- 
face centered at T disappears as well as the strong intra- 
band scattering responsible for the SDW ground stated 

To summarize, we have measured the spectra of 



Ba(Fci_ x NL r )2As2 as a function of doping and temper- 
ature in the SDW phase and found a broad distribu- 
tion of internal hypcrfinc fields consistent with long-range 
modulations of the local spin density around the dopant 
atoms in the ordered state. These observations are con- 
sistent with recent DFT calculations and neutron scatter- 
ing measurements! 14 i 17 The spin-lattice relaxation rate 
reveals strong 3D critical spin fluctuations that accom- 
pany the suppression of T/v with doping. Ni doping is 
qualitatively similar to Co, except that it contributes 
more electrons per dopant atomic We find no evi- 
dence for low frequency spin fluctuations in the ordered 
state associated with glassy behavior as observed in the 
cuprates i 16 ' 26 Future studies of the region of coexisting 
antifcrromagnetism and superconductivity may shed new 
light on the question of why the superconducting conden- 
sate can survive in the presence of pair breaking scatter- 
ing centers. 
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